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ABSTRACT 


The p (n, p)T reaction can be utilized in a proportional counter 
for neutron-energy determination in the region 100 kev to 1 Mev. This 
has already been demonstrated and the reaction cross section in this 
energy region has been measured with a counter. For practical 
application of the counter as a spectrometer, such as in health physics 
work, it is desirable to increase the counter efficiency. A proportional 
counter utilizing an anticoincidence ring to reduce wall effect has been 


constructed for this purpose. Some preliminary tests of the proposed 


system are described. 
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۶۶ئ) 
INTRODUCTION‏ 


lhe detection of neutrons and especially the measurement of 
their energy have from the first attempts presented a more difficult 
problem than for the other common nuclear particles. This is be- 
cause the usual methods depended on effects resulting from the charges 
of the particles. The experiments by Dee in which he investigated the 
ionization produced in air in a cloud chamber irradiated by fast 
ipee t were reported at the same time as Chadwick's announce- 
ment of the discovery of the neutron. Dee concluded that, if the 
neutron interacts with atomic electrons at all, the process produces 
not more than one ion pair per three meters of the neutron's path. 

If we accept an arbitrary description of neutron energy, E as 
intermediate neutrons: | kev to 500 kev, fast neutrons: 0.5 Mev to 
10 Mev, then we may say in general that the predominant reaction of 
intermediate neutrons with nuclei is elastic scattering, and that in the 
fast range many other reactions appear, the most important of which 
is inelastic — . Hydrogen- or methane-filled ionization 
chambers or proportional counters became an important method of 
detecting fast neutrons because the neutron can impart practically all 
of its energy to a proton in a head-on collision, and the recoil proton 
is then the particle causing the action of the counter. However, a 
determination of the neutron energy spectrum by the recoil method 
requires a thin radiator and a double collimation- -first of the neutron 
beam and then of the recoil protons--resulting often in too low an 
efficiency. Theoretically, the second collimation, may not be 
necessary because the neutron spectrum can be obtained from the 
recor! spectrum by differentiation, but this leads to large errors since 
we are effectively taking the differences of numbers of about the same 
size. Nuclear emulsions and organic scintillation counters have be- 
come very important detectors in neutron spectroscopy, especially for 


neutron energies above 1 Mev. 





For energies below lg@Mev it bee@ame apparent that essatacmimaic 
nuclear reactions might be useful for neutron spectroscopy. New 
methods were especially needed in health physics applications, where 
any schemes requiring collimation are automatically ruled out. Fast 
neutrons are of major concern to health physicists because of the 
difficulty of shielding against them, and their large relative biological 
effectiveness ^. In 1946 Feld proposed that a proportional counter 
utilizing the (n, p) reaction might be useful for neutron spectroscopy 
in the 10- to 1000-kev range, with the advantage over recoils of 
detecting lower energies and making collimation unnece Bam ni le 
considered the most promising reactions at that time: 


N14 (mee liés 


A ያቸ بت رو‎ 


However, asuhegexpected, thesegdo. not veryiwell fulfill the first of 
the general requirements, which might be listed as follows: 

(1) The reaction cross section should be fairly large and its 
variation with neutron energy smooth and free from resonances in 
the energy region of interest (and accurately known). 

(2) The Q value of the reaction should be small, because (E tQ) 
is the total energy measured. by a counter, and for accuracy E 
cannot be too small a fraction of the total. 

(3) There should be no low-lying excited states of the residual 
nucleus to cause an ambiguity in the energy of the reaction. particles, 

We can look for these requirements to be met over appreciable 
ranges only in light nuclei, in which the energy levels are more 
widely spaced, avoiding resonances and excited states. Also, for 
small Q values appreciable penetration of the Coulomb barrier for 
(n, p) or (n, a) charged-particle reactions would be limited to light 
nuelei. 

Two (n, a) possibilities are 


Li? (n, a)T, 
pom. a. 


=> 


The first has a large Q (4.78 Mev) and a resonance at E = 0.27 Mev. 
The second, widely used for thermal neutron detection, has an ex- 
cited state and fairly large Q values (2.79 and 2.31 Mev). 

With supplies of helium-3 becoming more available and a likely 
possibility being HESE p)T, helium-3-filled proportional counters 
came under consideration. In 1950 Coon showed that the cross 
section in the 1-Mev region was about 1 pargis , comparable to the 
scattering cross section of hydrogen. The thermal cross section for 
this reaction 1s large (5400 + 300 wh repay eee, These types of reaction 
cross sections obey the 1/v law for a wide range, but intermediate 
neutrons may depart from the 1/v law. 

Some work on hellium-3-filled proportional counters has been 
done by Batchelor and others at Harwell, وج‎ They used a 
small amount of helium-3 admixed with xenon or krypton as the 
counting gas and stopping gas to reduce wall effect. Also, a small 
amount of carbon dioxide was added to stabilize the multiplication 
process. In their final filling the pressure of krypton used was 164 
cm of mercury. Since they found that, for E = | Mev, a pressure 
of 8 atmos of krypton is required to reduce the wall effect to 10%, 
they set up a computer program to correct for wall effect. Their 
work remeasured the neutron cross section in the range 120 kev to 
1 Mev more accurately than previously, showing a much flatter 
energy dependence than the 1/v law. 

They also pointed out a basic source of ambiguity in energy 
determination other than counter defects. As we show in Chaper VIII, 
the maximum recoil energy of a He? nucleus is 3 O E The 
true energy counts run from Q for E -0to E, * Q. Thus, for 
3/4 E, = Q, the recoil spectrum begins to overlap the reaction 
spectrum. On the basis of Batchelor's value of Q as 770 kev, this 
means that ambiguities exist if neutrons of greater energy than 1.03 
Mev are present. The scattering cross section up to about 20 Mev 1s 
comparable in size to the reaction cross section of interest, but of 


course only a fraction of those neutrons scattered deliver a large 
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portion of their energy to the ضا‎ recols4lin any case, thespresence 
of high-energy neutrons reduces the accuracy. 

Gaseous counters have notoriously low efficiency, and as 
helium-3 has become more plentiful and less expensive, it seems 
worthwhile to investigate the feasibility of increasing the counting 
rate by increasing the quantity of helium-3 in the counter. As we go 
to higher pressures, special problems arise, such as the demands 
for gas purity, higher voltages, and better nase Also, the 
container must besmore massive, and proper consideration should 
be given to the effects of the masses upon the spectral distribution 
of the neutrons, primarily caused by inelastic collisions (trans- 
parency correction). For these reasons it would be desirable to 
use helium- 2 alone as the primary filling, utilizing it for the count- 
ing and stopping action. However, we soon find it very difficult to 
decrease wall effect to a reasonable value by helium-3 pressure 
alone because stopping powers decrease with lower atomic numbers. 
For instance, the following ranges of a 5.3- Mev alpha particle at one 
atmosphere have recently been measured 7 

He lima Al CAY 

krypton--3.0 cm, 

xenon--2.2 cm. 
For this reason it is desirable at the same time to determine the 
practicality of reducing the wall effect by surrounding the sensitive 
volume of the counter with an anticoincidence ring. 

Since helium-3 is quite valuable, a fairly elaborate filling 
system must be devised to insure no loss of gas in filling and sub- 
sequent recovery. The purity requirements are not so great as in 
noble gas scintillation counting, in which a small amount of contam- 
inant such as oxygen or other polyatomic gases may completely 
quench the useful scintillation property, but some purification arrange- 
ments will probably be needed in the filling system to achieve good 
energy resolution and to remove impurities which inevitably get in 


to the gas during filling and recovery. It is important that the 


Pa 





01 


[1 


loss of helium-3 be kept small during the purification process.  Purifi- 
cation procedures for noble gases using metallic calcium have been 
applied successfully for some time in removing all electro-negative 
impurities, nitrogen, hydPoPen, dnd tora We sser devPree car bep 6 
Also, the method of physical adsorption of impurities in an activated 
charcoal trap cooled to liquid nitrogen temperatures has been used 
successfully in a helium-3 system, and has the advantage of greater 


[11] 


problem due to the radioactive background from tritium contamination 


ease and convenience However, it appears there may be a special 
of helium-3, and this problem is increased by using higher helium-3 
pressure 

It was felt that useful information about these two techniques, (a) 
using helium-3 as the detecting medium and primary counting gas and 
(b) reducing wall effect with an anticoincidence ring, could be obtained 
by testing a counter with a filling of ordinary helium. The natural 
abundance of helium-3 is about 0.00013, so that it is not ordinarily 
noticeable. Owing to the large cross section at thermal energies, it 
can be detected in tank helium by use of a nuclear reactor as a source 
of neutrons, and this has been suggested as a way of surveying world 
helium supplies for helum-3 contri . Because it is a double- 
magic-number nucleus, me is very stable, and the detection of 
neutrons in a helium-filled counter is due to the recoil alpha particles. 
Helium-3 and helium-4 have the same electronic configuration, and 
we can obtain information about the counting characteristics of one 


by experimenting with the other. 
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CHAPTER II 


COUNTER SYSTEM 


A method of reducing wall effect in proportional counters is to 
detect the particles that escape from the sensitive volume and exclude 
them from the measurements. To accomplish this it is desirable that 
the walls of the main counter be transparent and completely surrounded 
by transparent-walled counters in anticoincidence. The feasibility of 
Geiger-Mueller counters with transparent walls -- so that wall thick- 
ness would not limit the radiation to that above a minimum-- was 
investigated in 1944, and it was found that symmetrically placed wires 


13] 


were constructed with good results. However, when two such counters 


could approximate a cylinder Single counters of triangular type 
were placed close together, there was a deleterious effect due to 
distortion of the field at the central wire. Better results were obtained 


with a regular hexagonal shape. 


COUNTER CONSTRUCTION 


A counter has been constructed in which a layer of gas surround- 
ing the main counter is in effect turned into a separate detecting device 
by a circular array of wires. The geometry is shown in Figs. 1 and ረ. 
The inner circle of wires, grounded to the cylinder, serves as the cathode 
for the main counter and as the inner cathode wall for the ring counter. 
In the outer circle of wires alternate wires serve as cathode and 
counting wires for the ring counter. Since there are 36 wires in each 
circle, we effectively have the main counter surrounded by 18 separate 
counters. 

An arrangement similar to this, referred to as a ''wall-less 
counter', has been applied to low-energy beta spectroscopy and proved 
to be a practical tue | lam, As pointed out, this prevents 
escaping high-energy particles from having an effect on the low-energy 
end of the spectrum, but the effective sensitive volume now depends 


on energy. However, inaccuracies because of this circumstance will 
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Fig. 2. "Gross Section ol counter wimlup, 





be much less than from uncorrected wall effect. Another application 
of this type of counter geometry is the so-called "Giles telescope" 
neutron spectrometer 115] „ In this arrangement the counter is long 
and narrow, and neutrons are collimated along the long axis. The 
requirement of a thin radiator is obviated by using gas recoils and the 
anticoincidence method to prevent registration of pulses that recoil 
oblique to the chamber axis. Monoenergetic neutrons then have a 
pulse-size distribution with a well-defined peak. A recent applica- 
tion of this scheme, utilizing hydrogen and propane fillings at 
various pressures, has resulted in energy spreads of 10% for neutron 
energies from a few hundred kev to 10 Mev 

Some of the constructional details of our counter follow. The 
counter cylinder has 1/16-inch-thick brass walls and is 4 inches in 
outside diameter. The total gas volume less the inlet piping is 1.93 
liters. The vessel has been pressure- and leak-tested with helium at 
20 atmospheres, making operating pressures up to 10 atmospheres 
possible within safety requirements. All wires are stainless steel of 
3-mil diameter. The central wire is supported by kovar insulators 
and is electrically available only at the base. Its length between 
thickened supports is 8.5 inches. If we assume that, owing to a fall- 
off in field strength near the ends, only 90% of this length defines the 
sensitive volume, the sensitive volume of the main counter is 888 cc. 
The ring-counter anode wires and all cathode wires are each strung 
from a continuous piece of stainless steel wire and threaded through 
holes in 3-inch-thick lava insulators at top and bottom, the anode 
wire system being electrically available through another kovar at the 
base. A major drawback to this method of stringing is the difficulty 
of repair in case of wire breakage. The volume defined by the ring 
counter is 600 cc. 

The anticoincıdence ring also cancels out particles that enter 
the counter from the outside or are emitted owing to contamination of 
the brass walls. Korff says a well-cleaned surface should emit 0.56 


[17] 


alpha particle per on per hour within a factor of ten and 
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Sharpe gives.0,09 alpha panmtilicle per a Ber hour for commercial 
a in normal operation the counter will probably not be 
sensitive to single secondary electrons knocked from the walls by 
gamma rays, but unwanted counts due to pile-up will be more likely 
in the ring counter, xheregphe density ofisecomea5my electzonsas val! bė 
greater. A proportional counter with walls of copper tubing 1/32 
inch thick has been used for counting gamma rays of energy up to 1.6 
0 The gamma-ray attenuation 'tenfolding length'' ıs 
approximately 60 کر‎ for all materıals in the energy region where 
Compton scattering is the most important process [20] This extends 
from less than one to several Mev for all elements. Multiplying by 
density, we find our brass wall is 1.34 e thick. The tenfolding 
length converted to a 1/e length is 26.1 سید‎ We have 
exp (-1.34/26.1) = 95%, therefore about 5% of all gamma rays imping- 
ing on the counter wall will be converted to secondary electrons. Of 
course, many of these electrons will be absorbed in the wall before 
reaching the counter gas, but gamma-ray counting efficiency would 
be reasonably high if the counter were operated at high gas 
multiplication. 

The voltage at which a counter under given conditions is operated 


determines the gas multiplication. For a given counter Rossi and 


Staub have shown the functional relationship of to ERI to be 





M = M 


where 
a = wire radius, 
b 


ዳን ም ول ات‎ ጋጋ 111.7 


1 


cathode radıus, 


V =voltage across counter. 
Since multiplication usually takes place in a cylindrical region of 
small diameter near the wire, we might assume our two counter 


systems identical for purposes of figuring multiplication, except for 


Ju 
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the difference in In b/a. The multiwire potential problem would be 
very dıfficult to compute, but if we take an average b value for one 
of the ring-counter wires (averaged over measurements to the inner 
circle of wires, to an adjacent wire in same ring, or to the wall), 
for the ring counter we obtain b/a - 225 and for main counter bla = 1000. 
The ratio of the logarıthms is 1.27, therefore for the same value 

of gas multiplication, thıs indicates that the main counter should have 
a higher voltage by 21 9ل‎ 127. Of course he determining 
factor on multiplication of the ring counter will be to make it as sensi- 
tive as possible in order to detect all particles that leave the main 
counter, but not so sensitive that a large number of background counts 


keep the desired counts from the main counter gated out. 


E BECTRONIC EIRGUITRY 


A block diagram of the electronic cırcuitry used with the 
counter is givenin Fig. 3. The blocks below the dotted line complete 
the neutron-spectrometer system, but investigation of the operation 
of the system was made without this by scaling both coincidences and 
anticoincidences. A list of units actually used is given in Table I. 
One channel of the dual-channel variable gate and delay unit can be 
used.for the variable- gate block, and the use of this unit for the 
pulse-shaping function could be replaced by simpler electronic cir- 
cuitry. This unit is very useful, however, for experimenting with 
different values of discriminator setting and of delay and length of 
output pulses. If two pulses arrive at the pulse-shaping block at 
about the same time, we want to insure that there is no anticoincidence 
pulse to gate the pulse-height analyzer on. To do this we delay 
channel No. 1 slightly so that the channel No. 2 output pulse will 
definitely be at the coincidence unit before the channel No. 1 pulse 
arrives.. Also, channel No. 2 pulse is made longer so that the 
coincidence unit will see it during the entire time that it sees the 


channel No. 1 pulse. The count on Scaler 1 * 2 indicates the number 
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Fig. 3. Block diagram of counter electronic circuitry. 
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# 1 


lables. . Units used in countempelieetr ani ¡cinc a 


Unit 
Northern Scientific Co., ser. 
UCRL 3 kv Divider Panel 
UCRL Mod 5, Dwg. IX5334 
UCRL Dwg. IX4353D 
UCRL Mod. 3, Dwg. 2T4404 
UCRL መ ህክ 3 


UCRL Variable Gate and Delay Unit 


Dwg. 2T8084F and 2T8154 


UCRL 8 Channel Quadruple Mixer 


Mod. 1 


Dwg. 3T5894-1 


UCRL 1024 Scaler, Mod. 2 
Dwg. 3T8934 


UCRL 2 Channel 1024 Scaler, 
Dwg. 4T6895 


alps 


Function 
Regulated High Voltage Supply 
Meter and Divider Panel 
Linear Amplifier, ch. # 1 
Preamplifier, ch. #1 
عم زر[‎ 310000777160۳ #2 
Preamplifier, ch. #2 
Pulse Shaping 


Coincidence Unit 
Scaler ۶ 


Scaler 1, ۸ 


د اوه مس مم æ‏ 








of pulses from channel No. 1 not registered by the pulse-height 
analyzer. The versatile variable gate and delay unit, used with one 
input pulse and in conjunction with the coincidence unit, can also be 


used as a single-channel pulse-height analyzer if desired. 


IONIC MOTION 


With this size of counter one may have to worry about the time 
lag between formation of the pulse by the main counter and by the ring 
Counter when ionization occurs close to'thne Sdgce"of the Main counter. 
We can compute the travel time for ions if we know the mobility. and 
how drift velocity depends on mobility and the field. Sharpe has 
developed convenient equations for computing ionic motion in 
TEE el In general, in cylindrical or spherical fields, the field 


18 





where 

y =b/a 

r = radius, 

V= voltage across counter, 

and for a cylinder n - 1. 
An ion crosses the system from one electrode to the other in a time 
given by 

Ts [a PH (yPP*1.1)] / [ps1] K (V/1n ሃ)"] . 
where p = ل1‎ for electrons and 1 for ions in the equation relating drift 
velocity w to mobility K, 
w = KeP 

We are concerned with the time for an electron to travel across the 
counter because most of the pulse 1s formed near the wire owing to the 
concentration of the field at the wire whether there is multiplication or 


operation as an ion chamber. If we take a case of V = 1000 volts, 


P = 2.5 atmospheres of helium, and use Sharpe's value of 
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Ke رو و‎ ioe مس‎ 21 - sec-Natmos , we find for electrons 


T = 7.55 usec. However, we expect that the mobility of a helium- 
carbon dioxide mixture may easily be greater by a factor of ten, 
based on the way CO, increases the drift velocity of electrons in other 
noble gases. This is explained for argon on the basis of the Ramsauer 
effect and the high first excitation potential for argon, 11.5 volts, in 
comparison with the low excitation levels in a CO, molecule [10] The 
first excitation potential of helium is even higher, 19.77 volts, 
Extensive work has been done on collection times in argon, krypton, 
and xenon, and it is reported that once a fraction of one percent of 
CO, is added, the effect of original polyatomic impurities in a noble 
gas may be assumed insigničicantLŽ2! „ The transit time computed 
above ıs reduced by the same factor as that by which the mobility is 
increased, and should not be too long for the coincidence requirements. 
Ihe transit time for fa Hes ion in pure helium to cross the system 
from wire to wall, computed for the same conditions as for the electron 
and on the basis of a mobility of 10.8 cm. EEE. 151850 sec. 


This time is not too ımportant, since the pulse can be shortened by 


shortening the amplifier time constant. 


E 


CHAPTER III 


ENERGY RESOLUTION 


A discussion of energy resolution depends on the complete theory 
of proportional-counter action. This has been well covered in several 


books and الى م‎ 10, 17, 18, 21, 23, 24, 25, 26] 


The attempt here is 
to point out important factors concerning the energy resolution of this 
particular counter, using developments from any of the basic sources 


but modifying the terminology as necessary in order to be consistent. 


INHERENT LIMITATIONS 


Inherent fluctuations in a proportional counter are from (a) the 
number of ion pairs released by monokinetic radiation, and (b) the 
size of the avalanche each ionization electron produces. From 
theoretical work based on the assumption that the first was a Poisson 
distribution (it was found experimentally to be smaller than this) and 
expanded upon experimentally, it was shown that relative variance of 


total size of avalanche is approximately given by 
E 
E: III- 1 


where 

V 7 variance (mean square deviation) of size of avalanche, 

mg = mean of number of ions initially liberated, 

m, = mean of number of ions produced in avalanche by one 

6000 

We note that this inherent spread is no longer a significant 
factor at the energies at which we will be working because in the 1- Mev 
energy range and at an assumed mean energy to produce an ion pair 
of 30 ev, mo is greater than 10*. These theories were useful, and 
experiments were done down in the kev. region where mg Was a smaller 
number. For our magnitude of Mo the relative variance would be 
around 0.01%. Since experience has shown that the relative spread of 
pulse height (relative standard deviation of pulse height from a mono- 


kinetic radiation) is at least 1% or 2% (owing to mechanical defects of 
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a counter or defects in gas filling) at these energies, the inherent 
spread is not a significant factor. The mean energy per ion pair, 

which determines mg, is expected to be practically constant over the 
(E, * O) range. This may not be true in very pure helium, in which 

the mean energy has been measured as 41.3 ev, in contrast to 29.7 ev 
when a small fraction of another gas is present. A possible explana- 
tion on the basis of charge exchange and the metastable states of helium 


1s given by Sharpe with reference to experimental اا ا‎ 


NEGATIVE-ION FORMATION 


An important gas defect is negative-ion formation due to electron 
attachment to a gas molecule. It does not occur in pure noble gases 
nor in N,, H,, 
the collisions of electrons with the halogens, بر‎ or water vapor. In 


CO, CO», or CH}. However, it is likely to occur in 


the development of Eq. (III- 1) the following equation was derived for 
the case where a fraction h! of ionization electrons reaches the wire 


without forming negative ions: 


i V | 
np mg - 
where "0 = variance of number of ions initially liberated. For h'= į 


(no negative-ion formation), this equation leads to the previous 
approximation that relative variance is about 1/mg when we use an 


experimental result that Vg ii pn d That is, the relative, variance is 


1( | هو 4 الا 
m, |3 hî ,‏ 


(A bibliographical review on the development of Eqs. (III-1) and (III-2) 


approximately 





is given in Reference 26). 
Thus, for accuracy in low-energy work it is recommended that 
the gas be purified to remove negative-ion-forming impurities. How- 


ever, since the factor in brackets in the expression for relative 


- 2 





varianceis multidied by the reciprocal of m,, we expect that in our 


case the purification requirements may e so stringent. This can 
be roughly checked by an approximate 5222801416) for the con- 
centration C of electronegative gas in a mixture in terms of the gas 
pressure P (in atmospheres), and Shubweg 6 (in cm) based on gas- 


discharge theory: 


-ረ 
وود‎ / ۰ 0 
Ce ከከከ (III- 3) 


The Shubweg 6 is the product of the mean free life of the electron and 
its drift velocity under the applied field. Values represenitative of 
the case in many gases were assumed for the electron mobility and 
probability of attachment per collision, h, for the gas. 

Suppose we wanted a gas purity that would make h' in our 
expression for relative variance equal 0.99. If an electron started a 
distance r from the central wire of a counter, the probability of its 
reaching the wire without attachment, by the definition of the 
Shubweg 6, would be =, Therefore we set 9 =,0499 04 
6 * 100r. Fora counter in which 3 cm would be a representative 
value for r and at a pressure of 2.5 atoms we get, from Eq. (III- 3), 

C = 0.01%, a fairly stringent requirement. If we try h' = 0.9 and 
make the same assumptions, we get C = 0.1%, or 99.9 purity 
required. Now, going back to the approximate expression for relative 
variance and assuming mg = 10”, we can solve for a value of h! that 
would give a relative variance of, say, 0.1%. This turns out to be a 
seemingly very easy criterion, h’ = 0.179. However, making the 
same assumptions as before, we find C is about 1%. But in the 
equation for C the assumptions for electron mobility and h, which 
exhibits resonances for some gases dependent on electron energy, 
could both conceivably be off by almost a factor of ten, and considering 
also other approximations, we could easily be back to our 99.99% 
purity requirement. The parameters in the gas-discharge theory are 
not accurately known for all the pure gases and cannot be predicted 
for mixtures, especially if the percentages are not very accurately 


known. 


وت 





RADIOACTIVE BACKGROUND 


Lhespr oblemi is. diffegentsifsthezeges ansanaaueciabi et ida 
impurity inthe helium-3 supply. The loss of energy resolution is then 
due to the pile-up of low-energy background counts due to the radio- 
active decay of the tritium. This decay has an end-point beta energy 
of about 18.9 kev [27] 

lo estimate the required purity of helium-3 with respect to 


21] 


Consıder the piling up of square pulses of uniform height W and equal 


tritium we must first analyze the effect of pile-up in general 


width 7T. A pulse height n W is produced if n pulses occur within a 
time T. lfa scaler discriminatpr level were set at n W, an erroneous 
count would then be registered. We can get approximate values for a 
counter by taking 7 equal to the equipment resolving time. 

If the occurrence mate ofi vents in the counter each giving rise 
to an energy E has a mean rate of ng: the mean number of events 
during a resolving time is noT and the standard deviation is ኣ 867 for 
ngT > l. Thus, the energy dissipated per resolving time is 
noT E NngTE. Now suppose the threshold energy of the system is 
such that it takes an energy of nTE to produce a spurious count. Then 
the following expression gives the number of standard deviations by 


which the number of counts during a resolution time must exceed the 


mean number in order for a spurious count to be produced: 
(nr - nor و۷8‎ = (m - ng) / Nag] NT. (111-4) 


From this we can estimate the probability Pg of receiving a spurious 
count during the resolving time for a given ng and 7. The spurious 
counting rate is given by I. = pisi. If we estimate T as 0.5 psec and 
wish to keep 1. less than 0.1 count per minute, p, must be less than 
I.T - 5/6- nmm 


the mean between 5 and 6 standard deviations. Let us require that 


This probability corresponds to a deviation above 


ng be less than n by 6 standard deviations. The average energy of the 


tritium beta spectrum is about 5.7 kev, therefore let E equal this 


- 10m 
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value. The range of beta particles of this energy is such that most 
of the beta energy will be spent in the counter gas, as shown in 
Chapter VI. Since the Q of the Des (n, p) T reaction 1s 770 kev, a 
value of nTE less than this will not be detected as a spurious count. 
Of course, a value of nTE that is an appreciable fraction of 770 kev 
destroys the accuracy of neutron-energy determination if enough of 
the pile-ups occur 1n coincidence with valid neutron counts. However, 
the number of accidental coincidences between two counting rates is 
proportional to the product of the rates and the resolving time; in 
this case the true and spurious rates are expected to be small enough 
that the problem of modification of true pulse-height counts will not 
be as serious as the number of spurious counts recorded compared to 


the number of true counts. If we set nrE = 770 kev, nis approximately 


Be p Using Eq. (III- 4), we obtain 
8 
a 3: 10 - np Jig. . 19762 | 
1/2 
"0 


Solving, we have ng 三 1.86 ' 10° disintegrations per second. 

If this is the allowable mean disintegration rate, we can solve 
for the allowable number of tritium atoms in the counter from the 
disintegration equation: 

aN dt =. (0.693/ T, />) N= - ng, 
where 
N = number of tritium atoms 
Ve = tritium half life = 12.26 years. 


This gives N = 1.03 ° ٧ tritium atoms, Or if each molecule is 


10 16 molecules. The 


diatomic, the allowable number is 5.65 ۰ 
sensitive volume of the central counter is Va = SCC, Adat ፉ 2.5 
atmospheres, the total number ofrmolecules is D P سس‎ = 5.66 ۰ ND 
molecules, where D - Loschmidt's number at 157 G mr 2:56 * 10! ?. 


This gives a tritium composition of about 0.0001%. This looks like a 
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very serıous purification problem if there is significant tritium contam- 


ınation in the helium-3 supply. 


SAT UR A DION 095 


It is apparent that for high-energy particles and high enough 
multiplication, saturation effects near the wire should eventually 
destroy the proportionality of the multiplication process. Experiment- 
ally, the ratio of pulse sizes has been found accurate to 1.5% for M 
less than a critical value of M_(E) at which saturation effects begin 281. 


Fr the counter used the experimenters found 
E - M_(E) = constant = 10° ev. )111-5( 


Using this for our counter application, where the maximum (E, * Q) is 


1.77 Mev, and allowing a safety margin of a factor of two, we have 
M, (E = 1.77 Mev) = 28. 


If we also check expected charge-to-length ratio in avalanche compared 
to charge-to-length ratio on the wire, a ratio of 4% corresponding to 
Eq (III-5), it appears that adoption of the equation for this counter 
should ensure no significant loss of proportionality due to saturation. 
High values of multiplication are also not desired, owing to increased 
voltage-stability requirements, because the curve of multiplication 
versus counter voltage in most cases starts out rather flat but soon 


rises into a steep ascent. 


POSITION OF IONIZATION 


Because of the long travel time of positive ions to the cathode, 
amplifier time constants must be used that ignore part of their con- 
tribution to the pulse. This makes the size of the pulse dependent on 
the position of ionization. A general rule is to use M >10 to avoid 
this effect when counter radius is 100 times wire radius. Let us 
investigate the effect for our counter, for y = b/a = 1000, and for no 


multiplication. On the basis of the expressions for fraction of total 


0 





induced voltage due to electrons and positive ions, Staub shows that 
the fraction of the volume of the chamber in which the ionization can 
originate so that more than half of the pulse is caused by the fast 
electron collection is b/(b*a), since the contribution to the pulse is 
evenly divided when r = N(a - b) is the position of original ioniza- 

xu eL Or this is yÁy*1), which equals 1000/1001, which seems to 
indicate that multiplication might not be necessary for large y values. 
This assumed position corresponds to r = 3l.6a for our counter. As 
a further check we assumed the ionization took place at r = 3l6a, 
and calculated the fraction of the pulse due to electron motion 
(ionization in about 90% of the volume would give larger pulses), 
using the assumptions about ionic motion used in Chapter II and 

some other equations from Reference 18. This fraction was 0.835. 
Amplifier time constants in reality would be long enough to register 
some of the positive-ion motion, which would improve this fraction. 
However, starting from r = 3l6a under the same conditions and 
assumptions, we calculate the ionic motion to contribute only about 
0.002 to the total pulse in a collection time of 5 usec. Thus we con- 
clude some multiplication is necessary to remove effects due to 


position of ionization. 


LID EFE EOL SG 


With no preventive measures there would be an abnormally 
high electric field at the ends of a counter wire. The loss in 
resolution thereby can be partially erased by having the wire 
shielded or thickened near the ends to give a low field. The variation 
of M along a counter wire has been measured and found to be import- 
ant at a distance on the order of a counter radius from the thickened 


[29]. 


ends In a high-pressure proportional counter designed to be 
used for neutron spectroscopy it was found that the pulse-height 
distribution varied radically over the extent of the counter in a pre- 


وات ادا 


pared with the resolution at the center, was significant when the 


liminary design without field tubes The loss of resolution com- 


Ao. 


investigators moved a narrow pencil of thermal neutrons with which 
they were probing the counter 0.5 inch away from the center. The 
wire was about 8 inches long. The simplest way of overcoming end 
effects is to use a long counter--that is, long compared with the 
radius. Of course, this is not always practicable. Our counter, 

with wire length 21.6 cm and counter diameter 7.62 cm, is nota 

very long counter, and we might expect considerable loss of resolu- 
tion due to end effects. The end effects can be eliminated by mechan- 
ical methods, which add complexity to the counter. A divided counter 
was constructed in which a bead divided the wire into two unequal 


[30] 


lengths Subtracting one spectrum from the other eliminates 
end effects. Field tubes extending over a guard tube eliminate end 
effects when field-tube potential is adjusted to a proper value and 

field and guard tubes are correctly located. Field tubes were proposed 


[32 


in 1951131) and used the same year to investigate the low-energy 


spectrum of tritium down to approximately 200 ev. 


WALL EFFECTS 


Wall effects as outlined for a gas recoil chamber by Rossi and 
Staub? are pertinent for charged-particle reactions, since in both 
cases the loss in resolution results when the particle that causes the 
ionization in the gas does not expend its full range in the sensitive 
volume. Three cases are: (a) The reaction takes place in the 
sensitive volume but one of the charged particles produced crosses 
the end boundary; (bj The reaction takes place in the sensitive 
volume but the charged particle hits the lateral wall; (c) The 
reaction takes place outside the sensitive volume, but one of the 
charged particles enters the sensitive volume. The anticoincidence 
ring does not help in those instances in which the charged particle 
crosses the ends. Thus, it is desirable to have a small ratio of area 
of ends to lateral area, or-- what is the same thing--a.long counter, 
as in reducing end effect. Pulses reduced in size by the wall effect 


contribute a continuum of energies from zero to (E + O) to the 
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pulse-height distribution. Wıthout the use of an anticoincidence ring or 
a magnetic field to curve the paths of the charged particles, the wall 
effect can be reduced only by increased size of the counter and ın- 


creased gas pressure. 


OLEE R EEEE OCES 


Inaccuracies ın energy resolution can also come from variation 
in applied voltage, variation in position and diameter of the wire, and 
dust particles ou the wire. These are discussed in the' general 
references. Any alone could be very detrimental, but it appears 


that all can be avoided with careful techniques. 
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GAS FILLINGS 


The counter was first tested by filling from a bottle of argon 
plus 4% carbon dioxide. The fılling and evacuation system used had 
an oil diffusion pump and liquid nitrogen trap, but no special purifica- 
tion procedure or trap was used in the gas-filling line. Satisfactory 
operation was obtained in counting gamma rays at pressures of one, 
three, and five atmospheres. Since the secondary electrons have a 
continuous distribution of pulse heights, we would not expect a flat 
plateau in the curve of counting rate versus voltage in the proportional 
region, but we do desire a slope that does not place too exacting 
requirements on voltage regulation. At 5 atmos we found about 3% 
increase in counting rate per volt in the 2600-to-2700-volt range on 
the main counter. This was at a high value of M; the range of pulse 
heights available was several thousand times the minimum detectable 
pulse height. At 3 atmos in the same voltage range, the change in 
counting rate per volt was about 3%. Of course, the M was greater 
Owing to the lower pressure. The ring counter at 5 atmos had about 1% 
increase per volt in the 2100-to-2200-volt range. Its counting rate 
was higher owing to its proximity to the wall. 

The helium filling was from a bottle of 99.99% pure helium; the 
desired percentage of carbon dioxide was previously introduced through 
a separate filling line. Both filling lines were run through dry ice- 
acetone slush traps when co, was used or a liquid nitrogen trap if 
helium alone was used. (CO, would be frozen out by a liquid nitrogen 
trap.) The use of the traps did not seem to have a large effect. The 
CO, was found very necessary to avoid spurious discharges in the 
counter., 

The stabilizing effect of polyatomic gases mixed with noble gases 
seems to depend on their ability to absorb photons, in contrast to the 
small-photon capture cross section of the noble gases (many of these 


photons are produced by metastable atoms), and aiso possibly on their 


m 





ability to receive energy from metastable states that have high 
excitation energies. Colli has investigated ultraviolet photons in the 
decay of metastable argon atoms and found them to be responsible, 
through a photocathodic process, for the starting and maintenance of 
corona currentin cylindrical argon روہ‎ d. The photons were 
produced in the Townsend avalanche on the wire. Another effect of 
metastable atoms is the release of secondary electrons from metal 

| Ihe positive ions can also eject secondary electrons 
from metals, the requirement being that the energy of ionization be 
greater than twice the work function of the metal, which is about 5 ev 
for common cathode walls. Helium is expected to be more troublesome 
than other gases because of its high ionization potential, about 24.5 ev, 
and the high energies of its metastable atoms, about 20 ev. Evidence 
for the difficulty of interpreting experiments with helium has been 


EE 


presented on the basis of the existence of the helium molecule-ion 
and metastable molecule. It had been shown the He, ion is present 
at very low pressures, and the concentration increases with pressure 
relative to the concentration of He 

The stabilizing effect of co, was demonstrated in the counter by 
determining the approximate maximum voltage on the main counter 
that could be applied without causing a rapıd breakdown when a radio- 
active source was placed a given distance from the counter. Total 
pressure was 5 atmos. The results are listed in Table ll. The break- 
down caused by background radiation alone occurs with pure helium 
at 800 volts. 

However, adding more CO, was not a solution because the 
pulse heights were thereby drastically reduced. It had been found 
previously that in mixture of argon plus 10% CO, there was a loss of 
pulse height above pressures of 3 atmos, probably due to columnar 
recombination, which is high in CO. This refers to recombination 
along the track of a densely ionizing particle. It was not possible to 


attain a large enough range of pulse heights above a minimum dis- 


criminator setting for counting recoils from a Pu-Be (a,n) neutron 
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Table Il; Percentage oi CO, i In versus voltage 
for breakdown on main counter. 


Percentage of CO Critical Voltage 





2 
4 1500 
6 2000 
1۳ 2500 
20 2800 
25 3000 
3361 3 — 


source without going to a counter voltage at which spurious discharges 
would occur. The situation improved with decreasing pressure. 

A final filling at 2.5 atmos with 2% CO, was chosen for further 
experimentation. The curve of coutitune ratesversissvoltage for count- 
ing neutron recoils from the Pu-Be source showed a slope of 0.84% 
change in counting rate per volt in the range of 1100 to 1200 volts on the 
main counter wire. For the outer wire there was a 4.3% change in 
counting rate per volt in the range from 800 to 1000 volts. The 
voltage-stability requirements at 5 atmos at a voltage at which the 
spurious discharges did not occur were not significantly different. 

One advantage of helium as the counting gas is readily apparent-- 
its reduced sensitivity to background gamma radiation. The range of 
tie secondary electrons in s/em“ is about the same, but helium has 
about 1/10 the density of argon. The low atomic number of helium is 
also an advantage in that the interaction of gamma rays directly with 


the gas is small. 
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سل للب 
COUNTING RATE AND EFFICIENCY‏ 


For anisotropic neutron flux in a proportional counter the count- 
- [2] 
g rate may be expressed as 
= mE TaD , 
Pe NV f (E) ቆሩ ጅ'ጫ) (E )dE, 
n 
where 
N = density of reacting nuclei Indc o 
iim volume of absorbing material m 
ሮ(ጅ,) = neutron cross section of absorbing nuclei 


Ly 


(cm رل‎ 
Pp = neutron flux per energy interval ا ےت‎ 
sec-erg), 
F(E) = detection efficiency defined as ratio of count- 
ing rate to rate of reacting events. 
Ideally, for an exoergic reaction and integral detection, F(E ) 13 
step function: zero for E « B and unity for E. > B, where B is 
threshold bias energy for equipment. This step function is modified 
by a geometric factor, G(E_), for wall effect as shown below. For 
differential detection, F(E) is ideally a difference between two step 


functions with thresholds B, and B,. For a recoil-detection process 


1 
in which, with no collimation, a neutron of any energy may give 

pulses from zero to a maximum, we must evaluate E Ee). which 
gives the fractional number of recoils giving pulses larger than Db. 


lhus, for a recoil detector, one has S 


7 “ء/3190 


The term FE) above is a correction factor for perturbation of 
flux caused by insertion of the absorber. This factor is neglected in 
the rest of the discussion. 

A solution of the equation for counting rate would in most cases 
be extremely complicated. Some insight can be gained by considering 
the simpler case of a mmoergic neutron flux. Then we can eliminate 


the differential and replace Pr with 中 (E ) in neutrons/cm^-sec. 
n 
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For neutrons of a given energy this leads to 


I= NV. o (Bt) ESPE (Ek). 


The efficiency of a counter is defined as I/q, whereq is the 
number of neutrons per second traversing the counter volume; 
q = o(E JA, where Ais the effective exposed area of the counter. Then 


comiter efficieney is 


N e c (E) ቆ(ጀጅ ) F(E) N Nr c (E) )م‎ ) 


=. L = 
የ ባጋ. መ Á 
If we call ١. عع‎ - d, the average distance traveled by a neutron in 
crossing the ern 1 counter efficlepssus 
where e, Noc (E. ) d is the efficiency of radiation defined as the 


R 
ratio of the rate of reacting events to rate of neutrons traversing the 


=No (E )dF(E |=e 


61377806-7 volume. 

The quantity N is readily calculated by expressing it as D P, 
where Dis Loschmidt's Number and P is the pressure in atmos. 

For the counter, we have d = 1.59 cm if we take the total sur- 
face area of the main counter as the effective exposed area for an 
isotropic flux. If we assume 2.5 atmos of helium- 3 and a neutron 
energy for which the (n, p) cross section is one barn, we have 
en” 0.01%. For the counter filled with ordinary helium at the same 
pressure and for a neutron energy for which the scattering cross 
section is 2 barns, we have ene 0.02% for recoils. 

The opposite extreme to an isotropic flux is a collimated neutron 


beam. Here the general equation for counting rate is 


DER f ነ: 87 Nuo) T, F(E )dE | 
n 


where 
A = area perpendicular to beam, 


T = thickness in cm of material parallel to beam. 


HG 





Assuming a mmoergic neutron beam as before, we can write 


| Peer (L 


[1 E ¿No (E) E 


and for No (E) ٣ل‎ p eR 7 No (E.) lo: 
In this case we see that the radiation efficiency of a cylindrical counter 
in a neutron beam is dependent on its orientation with respect to the 
beam. If the beam direction is perpendicular to the long axis of a 
counter of radius b, solving for the average thickness of path gives 

(1/2)b or 5.98 cm for the Main counter. Thus, the counter has a 
higher radiation efficiency under these conditions by the factor 


T/d = 3.76. Of course, e, would be even higher for a beam down 


R 
the long axis of a cylindrical counter. 


EXPERIMENTAL SETUP 


However, in experimental setups the counter usually is neither 
in an isotropic flux nor exactly in a beam. Consider placing a neutron 
source at point P, a distance c from the counter axis and opposite the 
center of its length, as shown in Fig. 4. If c is not too small and the 
counter dimensions not too large, it seems we could approximate the 
expected counting rate by assigning an average neutron flux to the 


counting volume: 


7 2 
ከርም 5 ۵۲ eff 


where Q = rate of emission of neutrons by the source and 


reff 7 N (c d), where d = bė 中 k 37 er . Then, 


I=NV 9 0/ d 
(assuming mE) = ] for the time being). The rate for neutrons 
traversing the counter in neutrons per second can then be expressed 


as this average flux times the area of the central rectangle, or 


- 1 
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Fig. 4. Geometrical arrangement of neutron source and 
counter. 
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= 2b: 2f = ds bf 
E 47 0 ۱ m ' cd 








This expression was set down to show that this method is equivalent 
to a common method of computing the solid angle subtended by a finite 
rectangular counter . For one of the small rectangles of area b f 


in the geometry shown the equation for solid angle given is 


tanı 2 215 1/0 
For a small angle we have tan 2* and for the entire rectangle the solid 
angle would be 

41 cad. 


Thus, we have 


Q. s(QU/ TT MIDI 7‏ )45 ۳ 4 د 
as above. Also one should note that we have assigned the same‏ 
radiation efficiency in this geometry as for a beam, as we can show‏ 
٧ ሀ V /4 bf (still neglecting F(E, ).‏ د 114 ع eg‏ 
f = (n/2)b, the average thickness‏ ظ 1/4 2 ۲۶۰ 2 ۶ ظ But we have Va a‏ 
found for a beam through the counter.‏ 

An arrangement like this was used with a Pu-Be neutron source 
of 62. = 1,54 + 10° neutrons/sec and d = 21.2 cm. Then we compute 
E.) = 1662 cpm for a pressure of 2.5 atmos of helium, and using 
an average scattering cross section of 2 barns corresponding to an 
average neutron energy of about 4 Mev. (Actually the Pu-Be neutron 
spectrum is complex from zero to about 11 Mev, and no attempt has 
been made to compute an accurate average value.) | 

If we obtain curves at different values of high voltage, as shown 
In Fig. 5, we can estimate the point Where they tend to zero pulse 
height. At this point د‎ ) is unity and the experimental value of 
counting rate should correspond to 1662 cpm. It appears to be only 
about 400 cpm, not a close correspondence, but there were so many 


inexactly known quantities and dubious averaging processes that it is 
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Fig. 5. Integral bias curves for main counter at various 
counter voltages, 
A. 1200 volts, 
B. 1165 volts, 
C. x35 welts, 
and using neutron source in arrangement of Fig. 4. 
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hard to draw definite conclusions. The scattering cross section is 
rapidly varying in this energy range and the manner in which the inte- 
gral bias curve approaches zero is not known. Also consider the count- 
ing rate of 175 cpm for the highest counter voltage used at ኳ pulse 
height of 10 volts, the lowest discriminator setting for which counts 
were recorded. If we let G(E, ) - 0.83--a value computed for an 
average alpha recoil in Chapter VI then, for a 10-volt bias, we have 
F(B/E,) - 0.127, which would indicate a detection of à small fraction 
of the recoil alpha particles. The curves for the ring counter, Fig. 6, 
seem to tend to a much higher zero-point counting rate, but itis 
believed that other events than the recoils are being counted, since the 


volume of the ring counter is less than that of the central counter. 


RECOIL-SPECTRUM CONSIDERATIONS 


It should be much easier to draw some definite conclusions about 
the suitability of the tube's multiplication factor for obtaining usable 
pulse heights if a monoergic neutron source were used, preferably one 
giving a recoil alpha in the energy range of interest, 770 kev to 1.80 
Mev. The counting results from a monoergic source would be much 
easier to analyze. 

This is illustrated as follows. For a monoergic source the 
neutron spectrum, N(E ) versus E is ideally a spike and N(E ) D 
the total number of neutrons; the differential pulse-height spectrum 
due to the recoils, n(W) versus W, where W is pulse height, is a 
rectangle; and the integral bias recoil spectrum, N(W) versus W, is 
ideally a constant-slope line. These are sketched in Fig. 7. These 
shapes of the differential and integral bias curves are justified as 
follows. The pulse height depends on the energy of the recoil nucleus; 
the counter and amplifier convert that energy into pulse height. The 
recoil nucleus, as shown in Chapter VIII, can have a maximum energy 
of K. E where K. = 4 2 di Dos ‚if we use| mass numbers for the 
masses. Let W be the energy of recoil nucleus (converted into volts 


of pulse height by the counter system). Then W can vary from 0 to 
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Fig. 6. Integral bias curves for ring counter at various 
counter voltages, 
A. 1025 volts, 
B. 975 volts, 
وت‎ 925 volts, 
and using neutron source in arrangement of Fig. 4. 
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Fig. 7. Ideal curves for monoenergetic neutron spectrum. 
A. Neutron spectrum. 
D. Differential bias recoil spectrum. 
C. Integral bias recoil spectrum. 
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K. E and the number of pulse heights from recoils from the monoergic 
source are spread out over this range of W values. The question is 
whether or not the distribution is uniform over the range, giving a 
rectangle. At the values of Ea 510061015066, And for eE, 
it is accurate to assume that neutron scattering is mainly S-wave and 
spherically symmetric in the center-of-mass system of coordinates. 
Wewish then to"express the*recoil'eWeTay ternas of E and the 
angle of scattering in center-of-mass coordinates (both energies in 
laboratory-system coordinates). Using conservation of energy and 
an eguation aeverogeal ?$- for relating the velocities of the neutrons 
in lab coordinates before and after scattering with 0, the angle of 
scattering in c.m. coordinates, we can show W = TE Ee (1-cos 0); 
then dW = 5 K, E, sin 0 dó. But in spherical coordinates, we have 
dð = 16/2۲ sin 6 , where dQ = element of solid angle, therefore 
CA E der (1/4). We can say Prob(W) dW = (0 (0)/0 لي‎ , 


where 
o (0) = differential scattering cross section (c. m. ), 
haw 8 integral scattering cross section, 
Prob( W) dW = probability for recoil nucleus to acquire an energy 
between Wand W+dW, 
(c (0)/c 92 - probability that neutron will be scattered through 
angle @ into dite 
It follows: Prob (W) = ٥ )6(/ ¢ : 4ፐ/ 35. E, E For an isotropic 
scattering distribution (0 ( ፀ)/ ። . = 1/4) and monoergic neutrons, 


Prob(W)is thus a constant, 1/K, E and n(W) is a rectangle of 
constant value proportional to N(E. ). The integral bias distribution 
is an integral of the rectangle, integrating from maximum W to zero, 
giving the triangular shape. 

For a non-momoergic neutron spectrum, N(E }--instead of being 
the total number of neutrons--represents the number of neutrons per 
energy interval, and we must integrate over the spectrum. We can 


say 


252 0 
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where K, 1s an arbitrary constant and وہ‎ (E 1s the elastic scattering 
cross section as à function of neutron energy. Thus. for a complex 
neutron spectrum the shapes of the differential and integral bias 
curves are very difficult to predict. For performing the integration 
in analytical form an expression for NIE, j and of E must be used 
Generally, for most sources NIE) goes from zeroat E 0108 
maximum and backto zero at E, - E (max). The intervening 
structure may be quite complicated. Let us assume 5 د‎ varıes as 


[18] 


for E, between 0.2 and 5 Mev, an approximate relation being Lo) F 


ت7 - 16 ٤‏ ۲ے 0 197۰5 Thıs ıs approxımately true‏ . ےا 
45/NE. for E, in Mev and E. in barns. itis less true for helium,‏ 
which has a peak of about 6.6 barns at 1.15 Mev and then tails off‏ 
Then if, for illustrative purposes, we assume‏ 


N(E )- E 3/2 sin E 
n n n 


where E. goes from zero to E, (many = tw Gn arbitrary units, we hame 


aspect shapeñthat can be readily integrated: 


E =r 
n 5 
n(W) = K J 1/(K E) م‎ E ^? ag 
1 a n n n n n 
E AE 
a 
= Ke (1 * COS (W/K y. 


pere W goes from 0 to K. T. To save writing W/K, as the variable 
٣٢ ነር next integration, let us change W from energy to arbitrary 
voltage units where W goes from 0 10 v. Since we have divided the 

W scale by K we must multiply the density distribution by K. The 
differential distribution 80.11. 


n(W) = K, (8 + cos W, 


= 38- 








The integral bias distribution is obtained by integrating this from 


the Tight: ۷ 
N(W) = - K, | (+ cos W) dW = K, (1 - W - sin W). 
TT 
The value of this distribution at W = 0 is the total number of scattering 
events, دی‎ This should equal the number of neutrons times the 


probability of scattering or 


Ki | N(E ) o (E) dE, = K, | E. sin E. dE, = TK. 
0 0 

The three distributions are sketched in Fig. 8. The sıgnificant 
trend is the way the integral bias distribution rises rapidly near zero 
pulse height, even with this assumed neutron spectrum weighted 
heavily toward the high-energy end of the spectrum. With a complex 
neutron-energy spectrum--as obtained, for example, from (a,n) 
sources --the shape of the integral bias curve and the manner in 
which it tends to the zero pulse-height counting rate becomes very 


difficult to predict. 
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Fig. 8. Assumed neutron spectrum and recoil spectrum 
for assumed conditions. 
A. Neutron spectrum. 
B. Differential bias recoil spectrum. 
C. Integral bias recoil spectrum. 
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CHAPTER WI 


سے 


PARTICLE RANGE AND WALL EFFECT 





Knowing the range of an alpha particle or proton in air, we can 
compute the range in helium if we know the relative atomic stopping 
power of helium compared with air. The range-energy curves to be 
used for alpha particles and protons in air at 15° C and 760 mm 


[38] 


pressure are those of Bethe, which, on the basis of experiment, 
corrected for low energies the original curves by Livingston and 

. L, The orıginal work was a theoretical treatment based on the 
Born approximation but also modified by experimental data then 
available. 

The basic equation giving energy loss per cm of path defines 

stopping power: 
REV AR (4x SN NB, 
where 

E = energy of incident particle, 

R = path length, 

e = electron charge, 

z = atomic number of incident particle, 

m = electron mass, 

v = velocity of incident particle, 

N = number of atoms per cc of stopping material, 

B = stopping number = Z In (2m 6 where 

Z = atomic number of stopping material, 

I -KZ-average excitation potential of atom of stopping 
material, 

K - Bloch constant (which is not actually a constant but varies 
somewhat for different atoms and is usually determined 
experimentally because it is difficult to calculate 
theoretically). 


The relative stopping power, compared with air, is defined: 


E - 





2 
一 = = ነ y 
5 B/Bg Z/Z yl (ln ረረ ln I „Vn 21۳7۷ - In Ly) d 


where the zero subscripts refer to air. Here s is seen to be a function 
of incident-particle velocity, and for different atoms increases more 
slowly than Z alone. We find it accurate enough for the purpose at 
hand to use this eguation, although many modifications to the theory 
can be made in computing B, as done by Hirschfelder and اہی ےون‎ 
for several different substances for protons. 

We take experimentally determined values of air from 
Reference 39: 16 2-3 5 ev and 7. O gern, ly = 44 ev 
m 0.328. To convert 


He 
the range in cm in air at NTP (15° C, 760 mm) to that in helium under 


[41] 
0 ۱ 
The calculation shows, for 1-Mev protons, s 


the same conditions we note that stopping power is proportional to 
N and B, therefore the range in cm must be inversely proportional to 


these quantities. The range in helium is 

R=R,' N/N . ጋጅ ۱ 
where the zero subscripts again refer to air and Sc s l. Now, N/N 
is approximately 2, since air is mostly diatomic. The calculation 


yields R = 6.13 Ro 
helium-4. For a l-Mev proton, Ro is 2.28 cm, therefore R is 14.0 cm. 


This should be fairly accurate for helium-3 or 


The range is inversely proportional to gas pressure, since density is 
proportional to pressure. Thus, at :2.5 atmos this proton range in 
helium becomes 5.60 cm, and at 5 atmos, 2.80 cm. 

To find the range of an ion of other isotopes of hydrogen (given 
the relation for protons) we note that of all terms in the stopping 
power formula only I is different for a given energy. Thus, a 
triton of energy E experiences the same stopping power as a proton 
of energy 1/3 E; but it goes three times as far since it has three times 
the energy. The maximum triton energy with which we will be con- 
cerned in the eS disintegration is about 1 Mev, and we figure its 
range to be 8.54 cm in 1 atmos of helium. Since this is less than a 
representative proton range from a Ee disintegration and represents 


an upper bound on the triton range, the order of magnitude of the 


=p = 


wall-effect difficulties can best be estimated by considering the 
proton range.  As'ointed out ንን ና. in explaining his computer 
program for wall-effect computation, the correlatıon between proton 
and trıton tracks plus a cylindrical geometry makes exact calculations 
impracticable; also the angular distribution of the ay Es p)T reaction 
1s not very well known. 

In integral bias counting involving the (n, p) reaction, detection 
efficiency is defined as the counting rate divided by rate of emission 
of protons. When the total range of particles is small compared with 
linear dimensions of a counter and distribution is isotropic, a value 


for detection efficiency often giwenms 
EME) 2 - Rp/2b 1 Rp y? e 


where we have an active cylindrical volume of radius b and length L, 
aud RE is that portion of total range of protons x ^x to Produce 
a pulse above the bias energy E of the discriminator „ The terms 
in the expression for F(E) reducing it from unity represent the 
magnitude of the wall effect and indicate the loss at the lateral bound- 
ary and end face, respectively. 

The expression for F(E) can be shown to be plausible as follows: 
Consider the fraction of the volume of a cylinder of radius b within a 
distance R of the side walls. This fraction is (2 bR - RIN 
Consider charged particles originating with equal probability along a 
radial line with ends defined by the outer volume. If any direction of 
emission is equally probable, one-half of the particles will be emitted 
with a radial component toward the center and not intercept the wall. 
For that half of the particles having a positive radial component, 
those originating on the inner edge of the fractional volume have 
practically zero probability of intercepting the walls (since only one 
direction of emission allows the particles to reach the wall), while 
those originating at the wall have 100% probability of reaching the wall. 
Assuming that this increase in probability is linear as we move out- 


ward toward the wall, we estimate that one-fourth of all the particles 


- UD 








originating in the fractional volume reach the wall. Then the wall 


effect is 1/4 time the fraction of the volume, or 


R/2b - m : 
For R << 2b, this is approximately R/2b. Similarly, the fractional 
volume within a distance R of the ends of a cylinder is 2R/L, and-- 
again taking 1/4 of this value-- we get R/2L as the wall effect due to 
toc ends: 

In our case the condition of large dimensions is not well ful- 
filled; also we assumed that the ionization from the particle other 
than a proton was confined to a very small region, which is not so 
true for a Ha reaction as for one involving heavier nuclei. In 
differential detection F(E) is ideally the difference between two step 
functions with different energy thresholds. "Therefore, the deviation 
from unity is a rough indication of the magnitude of wall effect on 
energy determination for a proton of given energy. For a 1-٧٠ 


proton in 2.5 atmos of helium, the geometric factor is 


۰۶۱ ۱۵۵ | RTS RL 
መመ with 2b = 7.62 cm, 2L = 38:9 cm, R = 5.00 cm, we have 
G(E)=0.12 or the wall effect is 88%. Doubling the pressure would 
reduce the wall effect to 44%, and at 10 atmospheres to 22%. 

For counting neutrons in ordinary helium we are concerned 
with the range of the alpha recoil. The maximum energy of the alpha 
particle is 16/25 E as shown in Chapter VIII. If we consider; the 
average energy of the Pu-Be neutron spectrum to be 4 Mev, this 
corresponds to a maximum alpha energy of 2.56 Mev, or an average 
E of about half this amount. Figuring the range of this alpha 
particle in 2.5 atmospheres of helium by our previous method, we 
get 1.38 cm. The wall effect is 22% for this range. Experimentally, 
a coincidence counting rate of about 16% of the total counting rate in 
the main counter was measured. This is a very rough check, but it 
at least appears that the ring counter is capable of reducing the wall 


effect. 


Jam 


ELECTRON RANGE 


To obtain the range agd nergy loss of electrons in estimating the 
effect of trıtium decays or background gamma counts, one can take the 
nearest tabulated value 8 of stopping power for helium as a function 
of electron eneTgy, given ip units ol Mgev por ee For at 
10 kev, we have 


 .‏ و وی ہد جس یا 


Likkuosvtake A xwes the average patheinebhe couwaber, 159 cm, stime suthe 
density of helium at 2.5 atmospheres, Ax = 7.53 - jor $ 06 
AE is greater than 16.9 kev (the product of - dE/dx ٠ Ax), because 
the stopping power increases with decreasing energy. For higher- 
enengy electrons it isplairly accurate toqwse"Feather"s"Rule.or 
Flammersfeld's الاك وو سوم‎ to find x and then multiply Ax/x times 


tHeveniesigysol theweleetron. 
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CHAPTER VII 
NEUTRON GENERATOR 


l'or testing purposes one desires a monoenergetic neutron 
source. The complex spectrum of (a, n) sources has been mentioned; 
the other class of neutron sources using radioactive nuclei are the 
photoneutron or (y, n) sources. These use deuterium or beryllium 
because other nuclei have thresholds above 6 Mev. To get high enough 
intensity the gamma source is surrounded so that there is an intrinsic 
energy spread due to different angles between neutron emission and 
incident gamma ray. This is small; for example, it is 3.37% for a 
Na T Be” 8.0 . However, this is not the main cause of 
energy spread. Owing to the large quantity of beryllium there is 
considerable neutron scattering in the source, and gamma-ray loss 
of energy by Compton scattering. 

It has been found possible to construct neutron generators 
utilizing the deuteron-deuteron reaction, D(d, fee > with useful 
neutron intensities without elaborate high-voltage equipment; for 
example, the generator of Zinn and Seely! 42] was built in 1937. 
Voltages on the order of 100 kv can give useful yields even though the 
reaction is not very efficient at this energy. With this low voltage we 
also have the advantage of almost monoenergetic and gamma-free 
neutrons. For a thick target the deuteron energy at the time of the 
reaction can vary from a maximum to zero. In the generator used 
in this experiment the maximum deuteron energy is 120 kev. Thus, 
using equations developed in Chapter VIII and Q = 3.28 Mev, we find 
that in the forward direction the neutron energy can vary from 2.46 
to 2.91 Mev. Ata laboratory-sybtem angle of 90 degrees the energy 
spread is from 2.46 to 2.49 Mev, or little more than 1%. The energies 
are weighted toward higher values because yield increases rapidly 
with deuteron energy. This neutron-energy range is convenient for 
checking the operation of the counter because the maximum recoil- 


alpha energy is about the same as the maximum (E t Q) energy we 
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wish to measure in helium-3, and about half the recoils are above the 


lowest sensitivity needed. 


PULSED- NEUTRON GENERATOR 


The neutron generator that is being put into operation provides a 
pulsed source of neutrons and is similar to the one described by 
Rubyl43! ExCCpDL Ior the method of mounting =the accelerater systempard 
the vacuum equipment. The accelerator system is described ın 
Reference 43; briefly, the deuterons are generated in a pulsed Phillips 
Ion Gauge type of ion source supplied with deuterium gas through a 
palladıum leak and accelerated to impinge upon a circular target of 
1-1/8-inch diameter. The targets are prepared by occluding deuterium 
in titanium. This type of target allows simpler equipment than one 
with heavy ice, by requiring no special cooling. It is planned to try a 
deuterated paraffin target to see if the yield is increased. The 
obuection o£ non-gaonoenergeticoneyutaonsaduesto the 29 8 re- 
action 1s not pertinent here because we are below the reaction 
threshold. The generator was designed so that adjustment of param- 
eters would not be critical; the yield is largely determined by the 
accelerating voltage and the condition of the target. (The reaction 
cross section decreases with deuteron energy, and the targets have a 
tendency to form a carbonous surface after a few hours! use. ) 

The pulsed operation is particularly convenient for experiments 
that require a reference time, such as in neutron-diffusion experiments; 
single pulsing has been used for cloud chamber experiments. In normal 
operation the ion source is pulsed at 60 cycles per second with 200- 
usec pulses. In this manner we can get an integrated yield, Q of 
more than 109 neutrons per secXrhüsuis gütyxcyclesf 00/2. 1955 
during a pulse the instantaneous neutron yield is Q /0.012 neutrons 
per second. For a cloud chamber experiment we used an external 
oscillator option of 400 cycles per second and increased the width of 
the pulse by an adjustment in the arc pulser to 500 psec. This is a 


duty cycle of 0.2. This mode of operation was gated on, when it was 
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desired to take a picture, for 100 milliseconds, giving a total pulse-on 
time of 20 msec per gate. Then the neutron yield per gate should be 
Q,/0r012 times 20 msec, or more than 10 neutrons in each burst, 
assuming the same effectiveness of the ion source when pulsed at 400 
cycles per sec. as at 60 cycles. It would be more realistic to calibrate 
at 400-cycles steady operation, but this increases the risk of damage 


to the ion source. 


NEUTRON YIELD 


In order to measure the integrated yield one should determine 
if the anisotropy of the yield is significant. It has been found that in 
the energy region below 0.5 Mev, the neutron yield as a function of 


angle in center-of-mass coordinates can be represented [44] by 


No) = Alo ON 
where Ais a constant and Bis a number that increases with deuteron 
energy. Thus, the yield is smallest at 6 = 90 degrees. The quantity 


B has been determined see down to 500 kev, and the 


[46] 5 


curve extrapolates fairly well into a theoretical determination 
low energies based on the D(d, p) Fr reaction. It is believed the 
angular asymmetry is the same for the companion reactions at 
corresponding voltages. We estimate Bas 0.6 at 120 kev. This 
means the neutron flux at Q = 0 should be (1 + B) times the neutron 
flux at 6 = 90 degrees. Or, performing an integration, we find that 
if Q. is determined by measurements at 9 = 0, we should divide this 
by (1 4 B/2+B)) for the true yield or 1.23 for B = 0.6. This is a high 
estimate of the anisotropy, because not all the deuterons have the 
maximum energy. With these small deuteron energies, we expect 
that the laboratory-system angles will not be much different from the 
center-of-mass angles. We need not outline the calculations here 

(a description of the two coordinate systems and some convenient 
equations may be found in Reference 47), but for 0 = 90 degrees. 

we find the lab angle 0 L“ 85.6 degrees for 120 kev, and for 0+ 0 
degrees, 0 = 94.4 degrees. 
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The generator yield has been monitored with a BF, proportional 
counter surrounded by 2.5 inches of lucite moderator and a thin sheet 
of cadmium. Care must be taken in its calibration since this 
arrangement does not give equal sensitivity over a wide range of 
neutron energies. For a particular counter we found a sensitivity of 
0.29 counts per neutron per m when a Pu-Be source was used, and 
0.37 counts per neutron per ea when a mock fission source with a 
smaller energy range of neutrons was used. However, it can be 
calibrated on the D-d neutrons by comparison with a calibrated 
Hanson and McKibben type el long counter. The lucite embedded 


counter is easier to place and adds less scattering objects to the 


experimental arrangement. 
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CHAPTER WIM 


NUCLEAR REACTION DYNAMICS 


3 


Consider the general case sketched above of a particle of mass 
M, and velocity V1 colliding with a target nucleus of mass M, assumed 
at rest. After the reaction the resulting particle of mass M3 flies off 
with velocity V3 at angle 8 The nucleus of mass M, has velocity 


v, at angle or. We will use nonrelativistic dynamics, mass numbers 


4 
for masses, and laboratory coordinates. 
First, write two equations for conservation of momentum for 


components parallel and perpendicular to ME 


MV, = M44 cos bi + مہہ‎ cos 01; (VIII- 1) 


M,v, sin 6, = M,v, sind, . (VIE 2) 


By transposing the last term of Eq. (1) to the other side, 
Squaring both equations and adding, we can eliminate 9, . Then 
replacing velocities by use of E = ž M wh for energy and rearranging, 


we obtain 


We can write another equation for E, {rom the definition of the 


Q of the reaction. Having assumed E, = 0, we write 


E 


E,=QtE,-E,. 


For the BE (n, p)T reaction, using Q = 0.77 Mev and similar 
symbolism, we write two equations for the triton energy: 


E en را ےر تہ‎ ሰሎ 0 
0 n p 


1 3 3 p 3 BZ 


AO + E E 
0 0 p 
With these two equations we can solve for the proton and 
triton energies for given values of E, and 91. 
For E = 0, the proton receives 75% of the O energy and the 


triton receives the remaining 25%. That is, 


E 
P 


Ed 


0.577 Mev, 
0.193 Mev. 


For E, = 1 Mev, a tabulation 1s given below. 


91 Triton 





(degrees) Ex % of 
(Mev) (E. FO) 

0 1.74 903 0.03 I 

45 L 87 656 025 14.2 

90 1.8 6.0 0.69 39 0 

5و 1.00 43.5 0.77 135 

180 0.67 37.8 1.10 o2? 


It is often desirable to have a convenient expression for the 
energy of the ejected particle at 0 ian 0 and at 90 degrees as a function 
of the bombarding particle energy. For the forward direction the two 


equations become 


5-0 


86”) ለ = ٣ 7 


waves «96 0 





LE,‏ - ربط + ل ۶ اط 


A solution for E, yields 


= 


1 


2 2 E 
Aea tata r n GA ANS 
where 


=M,Q+E, (M, - Mj), 


2A, (M; Mj) *4 E, M, M, , 


ል.'= & + My)“ < 


3 3 


For the B (n, p) T reaction this gives 





. 3 5 1 12 © 13 
د يك‎ 73 ۵ + 2750 ۳ E- 1 ۰ 
For the D(d, n) > reaction, 

3 1 1 6Q 5 


For 91. = 90 degrees, the two equations are 


- 14) M 
p - MEER Lo E., 
د‎ pe ed 


፡  -:ی‎ 


E 


Solving, we obtain 


E, ها + کی‎ + [(M,-M,)/(M, ተ እፅ )] E, 


For the two reactions above, at 90 degrees, 
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3 1 
E = O t 5 E. ! 
A 1 
E. 27 Q + 7 E : 
For the case of elastic scattering we set Q = 0 and M, = M ] 
therefore the two equations become 
E = - E. + ی‎ ው (E E M2 cos 0 
4 M 1 ٧ 3 ٧٧ (۳ ۳۰۵0 | jew 
4 4 4 
E, = E, - E, 2 


We are interested 10 ihe eneroy given to the recoil nucteus, 
therefore we eliminate E, from the first eguation by using the second. 
Then we set 0t, - 180 degrees for maximum E, corresponding to a 


backward recoil. Solving for E (max), we have 


4M M, 
E, (max) = ————— E 


d 
(M, *M,) 


For neutron scattering in helium-3, 


1 0 


- 2 
3 (max) = 2 E, ; 


and for neutron scattering ın helium-4, 


E (max) = — E 
a 
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